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We show the feasibility of a tomographic reconstruction of 
Schrodinger cat states generated according to the scheme pro- 
posed by S. Song, CM. Caves and B. Yurke [Phys. Rev. A 41, 
5261 (1990)]. We present a technique that tolerates realistic 
values for quantum efficiency at photodetectors. The mea- 
surement can be achieved by a standard experimental setup. 

03.65.Bz, 42.50.-p 



Optical homodyne tomography is a powerful tool to 
measure the density operator of a quantum system, with 
the possibility of detecting purely quantum features of 
the radiation field. Such technique was suggested in the 
past to reconstruct the density matrix of superpositions 
of classically distinguishable quantum states ||,|| , tradi- 
tionally known as Schrodinger cats, which represent one 
of the most celebrated examples of highly non classical 
states. Some experiments have been performed to detect 
Schrodinger cat states in atomic systems [[|. For radia- 
tion, a scheme to detect cat states has been proposed by 
Song, Caves and Yurke 0, however with no feasibility 
study for a real experiment, concerning in particular the 
main problem of quantum efficiency of detectors, which 
washes out the fringes visibility. 

In this paper we analyze a concrete experimental setup 
for the scheme given in Ref. and propose to detect 
the output field by means of homodyne tomography. We 
will examine the effects of quantum efficiency and present 
a method to compensate them allowing a good recon- 
struction of the Wigner function of the Schrodinger cat, 
and recovering the visibility of the experiment, even with 
quantum efficiency rj^ = 0.3 at the readout photodetec- 
tor, and r\h — 0.8 at the homodyne detector. Our pro- 
posal to reconstruct the Wigner function allows to ap- 
preciate the whole anatomy of the cat rather than just 
hearing its mew from the oscillations in a single quadra- 
ture probability distribution. To our knowledge, this is 
the first method for detecting the density operator of 
Schrodinger cat states for free radiation, feasible with 
current technology. 

Let us first briefly review the experimental scheme for 
generating Schrodinger cat states proposed in Ref. ||] 
(similar setups were later proposed in Refs. ). The 
main idea consists in feeding two orthogonally polarized 
modes of radiation, called "signal" and "readout" , both 
initially in the vacuum state, into a parametric amplifier 
followed by a half-wave plate. The parametric amplifier 
generates a correlated state of the two modes and the 
half wave plate rotates the polarization directions by an 



angle 0. The global state of the two modes at the output 
of this setup is given by 
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where V{r) = exp[r(asan- 



describes the action of 



the parametric amplifier, T(8) — exp[8(asa ! R —a ! s aii)] de- 
scribes the polarization rotator and the coefficients Bj tTn 
are given by 
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The rotation angle and the gain parameter r are related 
by the back-action-evading condition as sin 29 = 

tanhr. The following step of the scheme consists in de- 
tecting the number of photons at the readout mode. As 
a consequence of this measurement, given n r photons de- 
tected at the readout, the signal mode is reduced to the 
state 



where [^J denotes the integer part of x, 7t(n r ) is the 
parity of n r and P(k) is the probability of detecting k 
readout photons with a perfect photodetector, namely 
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In the scheme of Ref. || , after detection of the readout 
mode, the signal mode enters a degenerate parametric 
amplifier with gain parameter r s , described by the evolu- 
tion operator S(r s ) — exp[ir s (a| — a^ 2 )], which increases 
the distance of the two components of the superposition 
in the complex plane without changing the oscillating 
behavior of the number probability distribution. The fi- 
nal state of the signal is then described by the following 
quadrature probability distribution (the quadrature op- 



erator is defined as x$ = ^(a^e^ 
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where the eigenvalue X s((j>) of is the field quadrature 
component at phase 0, H n denotes the Hermite poly- 
nomial, A = [cos 4>(e~ 2rs cosh 2r cos 4> + i sin 0)] _1 and 
a = 1 + tan 2 <j) e 2r "/ cosh2r. 

In this work we propose to detect the Schrodinger cat 
at the signal mode using optical homodyne tomography 

. This technique consists in measuring the quadratures 
of radiation at several phases by means of a homodyne 
detector. The matrix elements of the density operator of 
the state are then measured as follows 



(6) 



where JC Vh represents the kernel function given in Refs. 
HD which depends on the value of the quantum effi- 
ciency of the homodyne detector rjh, while the overbar 
denotes the average over the experimental data at differ- 
ent phases. The behavior of the kernel function sets the 
validity limits of the tomographic reconstruction. These 
depend on the particular representation chosen to specify 
the density operator. In this work we always reconstruct 
the density matrix in the photon number representation: 
this is possible for any value of the quantum efficiency 
i] h > 1/2 (Ip). The rj h = 1/2 bound for the overall 
homodyne quantum efficiency is not a severe limitation 
in a real experiment, since good homodyne detectors can 
achieve values of rjh between 0.85 and 0.94 [p"i"|| . 

Let us now consider the effect of non unit quantum 
efficiency r\d at the readout photodetector. According to 
the Mandel-Kelley-Kleiner formula fL2| |, a detector with 
non unit quantum efficiency is equivalent to a perfect 
photodetector preceded by a beam splitter with trans- 
missivity r\d- Then, one can see that when n r photons 
are detected at the readout, the signal mode is left in the 
following statistical mixture of Schrodinger cat states 
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is the conditional Schrodinger cat state at the signal 
mode (||) evolved by the degenerate parametric ampli- 
fier, and Pr) d (k) is the probability of detecting k readout 
photons with quantum efficiency rjd, i.e. the Bernoulli 
convolution of the probability (|]). In Fig. [I] we plot a 
Monte Carlo tomographic reconstruction of the Wigner 
function of the statistical mixture ([t]) corresponding to 
an experiment with r\d = 0.3, rjh — 0.8, r — r s — 0.4 and 
n r = 2. As expected, the effect of non unit quantum effi- 
ciency is to smooth the oscillations in the complex plane, 
which are the typical signature of quantum interference 
(notice that the theoretical Wigner function is practi- 
cally indistinguishable from the one plotted in Fig. ||). 



Therefore, the resulting state is more similar to a classi- 
cal mixture of coherent states rather than a Schrodinger 
cat. The degradation effects on the cat due to non unit % 
can be seen also in Fig. |^, where the number probability 
for the same simulation is plotted: the probability still 
exhibits a non-monotonic behavior, but the even terms 
no longer vanish. In Fig. ^ we report a simulation of 
the quadrature probability distribution at <fi — 0, which 
would be seen following the original proposal ||, with 
the corresponding theoretical curve. 
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FIG. 1. Monte Carlo tomographic reconstruction of the 
Wigner function of the state in Eq. (Q), with r\d = 0.3, 
r\h = 0.8, r = r s = 0.4 and n r = 2. Data are collected 
for 70 different homodyne phases and 4 • 10 5 simulated data 
are used for each phase. 
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FIG. 2. Monte Carlo tomographic reconstruction of the 
probability distribution for the same state and same simulated 
data of Fig. [I]. The simulated values with the corresponding 
statistical error bars are superimposed to the theoretical val- 
ues (solid line). 




FIG. 3. Monte Carlo tomographic simulation of the qua- 
drature probability at (f> = for the same state of Fig. |l]. 
The histogram contains 20000 simulated data, while the solid 
curve is the theoretical distribution. 
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We will now present a method to compensate these 
dramatic effects of realistic quantum efficiencies. The 
main idea consists in the inversion of formula ([?]) , which 
gives 

x E ( k t 3 ) (1 - -yps^PvAk + j) ■ (9) 

j=o V / ™ d 

Hence, a generic /c-th Schrodinger cat component of the 
signal mode can be reconstructed by measuring all the 
signal states corresponding to different readout numbers 
of photons (larger or equal to k), weighting each event 
according to Eq. (^|) . In this way we have the additional 
advantage of using all data with n r > k than just those 
with n r = k of the plain detection in Fig. [j]. Moreover, 
by processing the homodyne data according to Eq. (||) 
we can reconstruct the whole family of Schrodinger cats 
\tps,k) for different k's at the same time. Notice that this 
method can be used for the reconstruction of any set of 
states |^s,fe) conditioned by an inefficient photodetection. 
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In our case convergence of the series (^) below the thresh- 
old % = 1/2 is due to the additional decaying factor 
Pr) d {k + j). An analysis of convergence of the series of 
errors as m Ref. @ shows that there is no lower-bound 
for rjd if r < \ ln(2 + \f%) ~ 0.658. Notice that the series 
convergence is slower for increasing r, which corresponds 
to more excited (macroscopic) Schrodinger cats. This 
implies that the more macroscopic the cat is, the higher 
77d must be in order to have a good reconstruction. 




FIG. 5. Photon number probability of the cat \ips,2) using 
the reconstruction algorithm ([]), with the same parameters 
as in Fig. The simulated values with the corresponding sta- 
tistical error bars are superimposed to the theoretical values 
(solid line). 

In Fig. U we plot the number probability for the same 
parameters of Fig. |[ the simulated experimental results 
with corresponding error bars are superimposed to the 
theoretical value. As we can see, the tomographic recon- 
struction is very precise and the oscillations of the prob- 
ability are perfectly resolved. In Fig. || we finally report 
the quadrature probability distribution for = super- 
imposed to the theoretical curve. The visibility is totally 
recovered, in contrast to the result in Fig. [|, which would 
have been obtained according to the original proposal of 
Ref. §. 



FIG. 4. Monte Carlo tomographic reconstruction of the 
Wigner function of the state \ips,2) with i]d = 0.3, r\h = 0.8, 
r = r s = 0.4, using the reconstruction algorithm (^). The 
same number of simulated data as in Fig. hi is used. 

In Fig. | we show a tomographic reconstruction of 
the same Schrodinger cat component of Fig. [I], with the 
same values of the experimental parameters, but using 
the reconstruction procedure based on the inversion (|^). 
As we can see, all the oscillations in the Wigner func- 
tion are properly recovered, and the destructive effects 
of low quantum efficiencies arc defeated. Notice that un- 
like other compensation methods based on the inverse 
Bernoulli transformation [ fl3|| , where the convergence ra- 
dius of the procedure is ij > 1/2, the present method 
works also for very low values of the quantum efficiency. 
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FIG. 6. Quadrature probability distribution at zero phase 
of the cat state \ijjs,2), with the same parameters as in Fig. ^. 
The simulated values with the corresponding statistical error 
bars are superimposed to the theoretical curve (solid line). 
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Let us now discuss in more detail the experimental fea- 
sibility of the scheme. All the devices needed in the ex- 
periment are available with the current technology. The 
parametric amplification can be realized for example by 
an ordinary KPT crystal pumped with the second har- 
monic of a mode locked Nd:Yag pulsed laser working 
at 80 MHz with 7 ps pulses fjj. The major problem 
encountered to detect Schrodinger cats in conditional 
measurement schemes is the cat's notorious fragility to 
any kind of losses and inefficiencies. The novelty of the 
present proposal is that using the reconstruction method 
based on Eq. (||) low values of rjd can be tolerated, 
and hence ordinary linear avalanche photodiodes with 
r]d ^0.3 can be used. On the other hand, the tomographic 
apparatus needed to detect the Schrodinger cat at the sig- 
nal mode is based on homodyne detection, with the pos- 
sibility of using high-efficiency photodetectors, because 
single-photon resolution is no longer needed due to the 
amplification from the local oscillator (LO). Moreover, 
the LO comes from the same laser source of the classical 
pump of the OPA, in order to achieve time matching of 
modes. In addition, since there is no fluctuating phase 
in the whole optical setup, neither in the second har- 
monic generation stage nor in the homodyne detection, 
the LO is also perfectly phase matched with the pump. 
The resulting setup is very stable and can take measure- 
ments for tens of minutes at a rate of 10 8 data/sec at the 
readout photodetector. The tomographic reconstructions 
presented in this paper were obtained with 2.8 • 10 7 ex- 
perimental data. In these examples the probability of 
detecting less than two photons at the readout photode- 
tector is ~ 0.9967. Therefore, taking into account that 
only the fraction 3.3 • 10~ 3 of experimental data collected 
at the readout is useful for the Schrodinger cat recon- 
struction, we can easily see that the whole set of data 
can be collected in a few minutes. Notice that increasing 
k, which corresponds to a more excited cat, the number 
of useful data decreases, and the reconstruction becomes 
slower. For example, to reconstruct the cat component 
with k = 4 only 2 • 10 -5 of the experimental data are 
useful, for k = 5 we can use only the fraction 1.7- 10~ 6 of 
data, and so on. Moreover, to reconstruct more excited 
cat components, higher index density matrix elements 
are needed for the Wigner function, and the effect of sta- 
tistical errors from tomography becomes more dramatic 
p5| , with the consequence that more data are needed 
to reach a prescribed accuracy. For these reasons, the 
more excited the cat is, the longer the experiment and 
the more difficult the state is to detect. Finally, it has 
been suggested Jl6) that by lowering the gain r some (non 
tomographic) homodyne visibility could be detected for 
high rjh (> 0.9) and for rjd as low as 0.3, even without our 
reconstruction algorithm. However, by lowering the gain 
the data acquisition rate is greatly reduced, whereas our 
method works also with high gains. 

In conclusion, we have shown the feasibility of a tomo- 
graphic reconstruction of a Schrodinger cat in an exper- 
imental scheme which is practical in a laboratory using 



standard technology devices. The problem of low effi- 
ciencies at the single-photon resolving detector, which 
was regarded as the major obstacle for experiments of 
this kind ||, has been solved by the implementation of 
a suitable tomographic data processing. The whole den- 
sity matrix of the cat, and hence all its characteristics 
(such as the photon number probability, the quadrature 
distribution and the Wigner function), can be measured 
in this way, whereas the plain homodyne detection pro- 
posed in the original scheme of Ref. || would not have 
provided visible probability oscillations with the available 
low-efficiency single-photon-resolving detectors. 

We thank T.F. Arecchi for illuminating discussions on 
the experimental setup. This work has been supported 
by the PRA-CAT97 of the INFM. 
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